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Abstract 

The polarization plane of the cosmic microwave background radiation can be rotated either 
in a magnetized plasma or in the presence of a quintessential background with pseudoscalar 
coupling to electromagnetism. A unified treatment of these two phenomena is presented for 
cold and warm electron-ion plasmas at the pre-recombination epoch. The electron temper- 
ature is only relevant to the relativistic correction of the cold plasma results. The spectrum 
of plasma excitations is obtained from a generalized Appleton-Hartree equation, describing 
simultaneously the high-frequency propagation of electromagnetic waves in a magnetized 
plasma with a dynamical quintessence field. It is shown that these two effects are compara- 
ble for the plausible range of parameters allowed by present constraints. It is then argued 
that the generalized expressions derived in the present study may be relevant for direct 
searches of a possible rotation of the cosmic microwave background polarization. 



If linearly polarized radiation passes through a cold (or warm) plasma containing a mag- 
netic field, the polarization plane of the wave can be rotated since the two circular polar- 
izations (forming the linearly polarized beam) are travelling at different speeds. This effect 
has been studied in a variety of different (but related) frameworks even in relativistic QED 
plasmas and in the presence of extra-dimensions (see, for instance, [1]). 

The cosmic microwave background (CMB) has a weak degree of linear polarization, which 
is the direct result of Thompson scattering. If CMB is linearly polarized, then its polarization 
plane can be rotated provided a sufficiently strong magnetic field is present around the time 
of decoupling [2]. Typical values of magnetic fields, compatible with other astrophysical 
constraints, are of the order of B < 10~ 4 G at the decoupling time (see, for instance, [3]). 

An apparently unrelated possibility is that the quintessence field a has pseudoscalar 
couplings to electromagnetism [4]. Couplings to ordinary matter, even if suppressed by the 
Planck scale, may lead to observable long-range forces and time dependence of the constants 
of nature. An approximate global symmetry may suppress these couplings even further as 
argued in [5] (second reference). This possibility may also imply a rotation of the plane of 
polarized emission, and radio-astronomical implications of this type of cosmic birefringence 
were investigated through various steps [5] . Typical values of the mass of o are of the order 
of 10~ 33 eV in such a way that, between redshifts and 3, a can start dominating the 
background with energy density mVg ~ A 4 , where A ~ 1CT 3 eV and ctq — Mp ~ 10 18 GeV. 
The effect of parity violating interactions on CMB polarization has been recently discused 
in [6]. 

In this paper a unified discussion of Faraday rotation and cosmic birefringence will be 
presented for the specific case of the pre-recombination plasma. The possible interference of 
these two effects will be scrutinized in both the cold and warm plasma approximations. A 
related aspect of the present analysis will be to study the range of validity of the Faraday 
rotation estimates. 

Let us start by discussing the typical scales involved in the problem. Right before decou- 
pling, the temperature of the plasma is of the order of 0.3 eV and the typical free electron 
density can be estimated as 

n e ~i c (fi b ^)(l + z) 3 xl0- 5 cm" 3 , (1) 

where x e is the ionization fraction, /i is the indetermination associated with the Hubble 
constant, Qb is the fraction of critical density in baryons, and z is the redshift. For typical 
values of the parameters the electron density is of the order of n c ~ 10 3 cm~ 3 . 

Right before decoupling, the plasma is globally neutral and the ion density equals the 
electron density, i.e. n\ ~ n e = n , denoting with n the common electron-ion number 
density. The global neutrality of the plasma occurs for typical length scales L>A D where 

A ° = \Z3k " 10 ( W^) (blTev) V2 ^ (2) 
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is the Debye screening length whose value is much smaller, for the same choice of physical 
parameters, than the electron and photon mean free paths, i.e. £ c ~ 5.7 x 10 7 cm and 
£ 7 ~ 10 4 (1 + z)- 2 ^^)- 1 Mpc. 

If the plasma is not magnetized, the only relevant frequency scales of the problem are 
the plasma frequencies, which can be constructed from the electron and ion densities, i.e. 

/ 47re 2 n n ( n e \ Aite 2 n ( n e 

^ = V^T"Hl0^n^J ^ = V^^" 40 ll0^n^J ^ 

(3) 

where, for practical estimates, the numerical value m; ~ m p = 0.94 GeV has been assumed. 

The frequencies given in Eq. (3) enter the dispersion relations, determining the group 
velocity of an electromagnetic signal in the plasma. The plasma frequencies for both electrons 
and ions are much larger than the collision frequencies constructed from the inverse of the 
mean free paths, i.e. u cc ~ 520 Hz and uj ci ~ 2.4 Hz. The plasma can then be described, to 
a very good approximation, within a two- fluid framework [7, 8]. 

If the plasma is magnetized, two new frequency scales arise in the problem, namely the 
electron and ion gyrofrequencies, i.e. 

u, Be = ^ c 18.08 (-§-) kHz, u Bl = ^ * 9.66 (-§-) Hz, (4) 
m e c V 10 d G/ m ; c \10 6 G/ 

where Bq is the magnetic field strength at the corresponding epoch. 

The electron and ion gyrofrequencies, together with the plasma frequencies of Eq. (3), 
affect the dispersion relations in the case of a magnetized plasma. The maximum of the CMB 
can be determined, from the Wien law, around 1.7 GHz; we will then be interested in the 
frequency range between, say, 1 and 100 GHz. In this frequency range wcmb uj pc > ub c - 
Since u pe > u p \ and ub c > wbu we will also have wqmb ^ P i > ^Bi for the fiducial set of 
parameters introduced so far. 

Defining the appropriately rescaled electron and ion densities, n e = a 3 n c and n ; = a 3 n ; in 
a conformally flat geometry of Friedmann-Robert son- Walker background characterized by 
a scale factor airf) and by the line element ds 2 = a 2 (r])[dr] 2 — dx 2 }, the continuity equations 
for the charge densities reads 

n' c + 3w e Hn e + (w e + 1)V • (n c v c ) = 0, (5) 
n[ + 3w{Hni + (w, + 1)V • (mi*) = 0, (6) 

where Ti = a'/ a; the prime denotes a derivation with respect to the comformal time coordi- 
nate rj; w is the barotropic index for the electron or ion fluid. 

Both electrons and ions are non-relativistic before decoupling. Hence the barotropic 
index w will be close to zero to a good approximation. For instance, the energy and pressure 
densities of an ideal electronic gas are given by 

p e = n e (m e + ^T^j , p e = n e T e , (7) 
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and since w e ,i — T e ^/m e ^ w e ^ <C 1 as far as T eji <C m^ e . 

In the cold plasma approximation the temperature of the ions and of the electrons van- 
ishes. In the warm plasma approximation the temperature of the two charged species may 
be very small but non-vanishing. The warm plasma treatment will lead, in practice, only to 
an effective correction of the plasma frequency. Since the cold plasma results turn out to be, 
a posteriori, rather accurate, the discussion will be presented in terms of the cold plasma 
description, while the results of the warm plasma treatment will be stated later. 

To have a self-consistent set of two-fluid equations, Eqs. (5) and (6) will be supplemented 
by the evolution equations of the velocity fields and of the electromagnetic field, namely 

PeK + nv e + (v e ■ V)u e ] = -n e e(E I^xfl), (8) 
+ m + (vi ■ V)«i] =iiie(£ + ^x5), (9) 

— » _ — * — * _ — * 

where E = a £ and B = a 2 B are the conformally rescaled electromagnetic fields obeying the 
following set of generalized Maxwell equations: 

V • E = Aneim - n c ) + ^Vtx • B, (10) 

V-S = 0, VxE = --B' (11) 

c 



V x B = -E' - ^ 

c M 



a'B + Vtr x E 



Aire, _ . . 

H {riiVi- n e v e ), (12) 



where the coupling of the electromagnetic field to the quintesssence field a can be derived 
from the action 



' a = J d A x^-g^-g^d a ad p a - W(a) + -]^f<#F* 



(13) 



where F a p and F a @ are the Maxwell field strength and its dual; g a p is the metric tensor. 
In the limit a — ■> 0, Eqs. (10)-(12) are the usual curved-space two-fluid equations [3]. The 
evolution of the quintessence field will obey the homogeneous equation 

dW 

a + 3Ha + — = 0, (14) 
oa 

where the overdot denotes a derivation with respect to the cosmic time coordinate whose 
differential is related to the conformal time coordinate as dt = a(rj)dr); it will be assumed 
that 2W(a) = m 2 a 2 . From Eqs. (10)— (12), the presence of the quintessence field introduces 
a further frequency scale into the problem, namely uj a = c(f3/M)a' . If the quintessence field 
dominates today (or between redshifts and 3), <j(t ) ~ A 2 /m where, typically, A ~ 10~ 3 eV, 
and m ~ 10~ 33 eV. The value of a can be estimated, today, from Eq. (14) and it is 
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a(t ) ~ m 2 /H . Hence, recalling that prior to quintessential dominance, & scales as a 3 , 
from the previous expressions 



/H dcc \ 2 mA 2 o(M\ 6 



0"dc 



where H dcc ~ 10 6 H ; the values of A and m are the ones discussed above with M ~ M P . 
The value of (3 is rather uncertain and a conservative limit from radio-astronomical analyses 
would imply /3<10 — 3 for 0< z ^1 [5]. Notice that the action (13) can also be relevant in a 
class of baryogenesis models [9] (see also [10]). Similar actions also arise in the framework 
of theories with varying coupling constants (see for instance [11]). 

Consider the combined effect of the homogeneous quintessence field (i.e. Vcr = 0) and of 
a background magnetic field on the spectrum of plasma excitations. Equations (5) and (6) 
together with Equations (8) and (9) and (10)-(12) can then be linearized in the presence of 
the weak background magnetic field B , i.e. 

n e>i (r),x) =n + Sn eji (r],x), B(r),x) = B + 6B(r),x), 

v ,i(v,x) = 5v Cji (rj,x), E(r],x) = 5E(r],x). (16) 

Using Eq. (16), the system of Eqs. (5)-(9) and Eqs. (10)-(12) can be written as 

5n' c + n V • 5v c = 0, 5n[ + n V • 5v { = 0, (17) 

(18) 



6% + H6v e = -— 
m. 



c 



e 



5E + — x B 

c 



5v{ + HSvi 

rrii 

V x 5E = --SB', V ■SE = Aire{Sn i -Sn e ), (19) 

V x 5B = -5E' - 4-ytB + ^^(Svi - Sv e ). (20) 

c M e 

From Eqs. (17)— (19) the relevant dispersion relations and the associated refraction in- 
dices can be obtained by treating separately the motions parallel and perpendicular to the 
magnetic field direction. Defining the current direction parallel to the magnetic field as 

— # 

j\\ = noe(6v it \\ - <Ju e ,||), Eqs. (18) imply 

Jii+nl, = ^1 + ^1)6^, (2i) 

Since the variation of the geometry is slow in with respect to the typical frequencies of plasma 
oscillations, the following adiabatic expansions can be used: 

J\\(v,x) = nAfy~ ir *' uW) > <^llM = SE^e-'r^-W). (22) 
Thus, defining a = iH/uo <^ 1, Eq. (21) implies that 

J ^ = ^u(l + a) 5E ^ (23) 
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Inserting Eq. (23) into the parallel component of Eq. (20), the following equation can be 
obtained: 

(V x 8B U )\\ = -£e\\(u,a)8E\^ - ^'5B loj , (24) 
where the parallel dielectric constant is 

2 2 

u 2 (l + a) u 2 (l + a) 

Notice that only the leading curved-space correction has been kept. 

With a similar procedure, also the motion in the plane orthogonal to the magnetic field 
direction can be solved and the evolution equations of the electric and magnetic fluctuations 
can then be written, in compact notation, as 

Vx5E LU = i^SB, (26) 

Vx^ w = -£e(a, u)5E u - ^SB U , (27) 

— * — * 

where 5E U and SB^ have to be understood as column matrices containing, in each row, the 
components of the electric and magnetic fields in each of the three spatial directions, while 
e(co, a) is a 3 x 3 matrix given by 



e(co, a) 



I €i(u,a) ie 2 (ui,a) \ 

— ?e 2 (cj,o;) ei(u,a) 
V e\\(u,a) J 



(28) 



where €||(u;,a) is defined by Eq. (25); e li2 (^,a) are instead 

ei(w, a) = 1 1 — , (29) 

uj 2 (a + l) 1 - ^ (a + 1) - 



2 2 

e 2 {u,a) = — — \ — — — . (30) 

uj u 2 (a + I) 2 — cjg c to uj 2 {a + I) 2 — cu^ i 

The coordinate system can be fixed by setting k x = and k y = ksm9, k z = kcos9 with B 
oriented along the z direction. Since Eqs. (26) and (27) imply 



VxVx^, = ^e(cu, a)6B u - -^VV x 6B U , (31) 



the Fourier transform of Eq. (31), in the coordinate system selected previously, leads to the 
generalized Appleton-Hartree equation: 

/ [1-^1 -^ + ^(6)] 
^ 8B^ =\i[% + <£40)\ [c 2 (0) - m s{9)c{e) | | 8B k ^ y \ = 0, (32) 

V -i%8{0) -8(0)c(0) [s\9)-^\ 
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where the refraction index n = c/v has been introduced so as to eliminate the comoving 
momentum, in such a way that k — to/v — noo/c; we have written c(9) = cos 9 and s(9) = 
sin 9. From Eq. (32), A^ = A, where the dagger denotes the transposed and complex 
conjugate of a given matrix. 

The non-trivial solutions of the system of algebraic (homogeneous) equations given by 
formula (32) comes from setting the determinant of the coefficients equal to zero, i.e. det^4 = 
0. After some algebra it is found that the determinant vanishes if 



s 2 {0) 



-c 2 {9) 



1 



rr 



2 ej 



2Ve L 
1 1 

~ 2 ~^ 



,2 



2n 2 co 2 



1 1 
— + — 

£l £r 



n 



oo„ 



n 2 uo 2 e R e R 



+ 



CO n 



u J oo \€l 



1 



1 

£r 



c{9) = 0, (33) 



where the right-handed and left-handed dielectric constants have been defined as 



CR — Cl + ^2 



— — 62 



Pi 



'pe 



ou[u(a+ 1) — ub\\ oo [u(a + 1) + u Bi 



CO 



CO. 



(34) 

1 Z21 . (35) 

io[uj(a + 1) - lo Bc ] io[uj(a + 1) + to Bi ] ' 

If co a = 0, then Eq. (33) reduces exactly to the Appleton-Hartree equation known from 
two- fluid plasma theory [7] , with the minor difference that the leading dependence upon the 
background geometry appears in e Rt ^ through the function a . The dispersion relations for a 
wave propagating parallel and perpendicular to the magnetic field direction can be obtained 
by setting, respectively, 6 = and 9 = 7r/2 in Eq. (33). Consequently, the relevant equations 
determining the refraction index are, in this case 



(n 2 -e R )(n 2 -e L ) = 0, 
(n 2 -e||)[n 2 (e L + eR) -2e L e R ] = 0, 9 



0, 

TT 

2' 



(36) 
(37) 



Equation (36) gives the usual dispersion relations for the two circular polarizations of the 
electromagnetic wave, i.e. n 2 = e R and n 2 = e L , while Eq. (37) gives those for the "ordinary" 
(i.e. n 2 = e||) and "extraordinary" (i.e. n 2 = 2e R eL/(eR + plasma waves [7, 8]. 

In the presence of a dynamical quintessence field weakly coupled to electromagnetism, 
uo a 7^ and the resulting dispersion relations are, according to Eq. (33): 



n 



n 4 - 



00a 

— n - e R 

CO 



2 , 00 (j 

n H n 

CO 



fLfR . ( 00 (j 

h en + — 

L £l + e R V co 



n 2 + 2 



0, 

CllCLeR 



0, 



9 = 0, 



(38) 
(39) 



£r + £l 

From Eq. (38) the generalized Faraday rotation experienced by the linearly polarized 
CMB travelling parallel to the magnetic field direction can be obtained as 



oJ 

A$ = — 

2c 




AL, 



(40) 
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where u ~ wcmb at decoupling; AL is the distance travelled by the signal in the direction 
parallel to the magnetic field direction. It is interesting to compare the contribution of the 
terms depending upon uj g and those depending upon the background magnetic field intensity, 
i.e. the terms appearing in the squared brackets. 

According to Eqs. (3), (4) and (15) the leading contribution to the generalized Faraday 
rotation arises as the sum of two dimensionless ratios, i.e. 

p^wj^vja-N* (41) 

Using Eqs. (3), (4) and (15), it can be verified directly that the two contributions appearing 
in Eq. (41) are of the same order, provided 



io-y \m p ) ~ uo- 4 g; uo- 3 cm 3 



(42) 



where in the left-hand-side we computed the quintessential contribution, at the right-hand- 
side the contribution of the magnetic field has been reported. The values of (3 and M 
are, as previously stressed, conservative estimates, compatible with present bounds [5]. For 
the typical (allowed) parameters reported in Eq. (42), the two contributions then are, as 
anticipated, of the same order. For the explicit value of wcmb the numerical value wcmb ~ 
THz has been used (notice, in fact, that the maximum of the CMB is of the order of the GHz 
today and of the order of the THz around decoupling). Given the possibility of absence of 
hierarchy between the two dimensionless ratios, Taylor expansions of the obtained disperion 
relations should be treated with care, especially in view of possible experimental applications 
of these findings. In a complementary perspective, when analysing the possible rotation of 
the CMB polarization, it seems then preferable to adopt the generalized formulae derived in 
the present study. 

Concerning the obtained results, a few comments are in order. 

• the two contributions have a frequency dependence that may allow us to disentangle 
their relative weight since the magnetic contribution vanishes for u > uqmb as 1/cj 2 ; 

• for the set of parameters chosen in Eq. (42) the typical generalized Faraday rotation 
measurement (i.e. A$ = RM \q MB , A being the wavelength) will be of the order 
of 200 rad/m 2 if we use the fact that the optical depth for Thompson scattering at 
decoupling is / x e n e a(rj)dr] ~ a^ 1 where <tt is the Thompson cross section; 

• situations can be envisaged (f3 <C 10 -3 ) where the quintessential contribution is far 
smaller than the "magnetic" contribution; 

• the dispersion relations for radiation propagating orthogonal to the magnetic field 
intensity also lead, in the case oj a ^ 0, to a rotation where the magnetized contribution 
is absent and A<3> cx w ff AL; 



7 



• in the (generalized) extraordinary dispersion relations (see Eq. (39) a resonance may 
appear, when, (3 <C 1CT 3 and for uo r ~ \J^p C + ^ic! however uj v <C wcmb; 

The derivation presented so far refers to the case of cold magnet oactive plasmas. The 
effects arising from the finite temperature of the electrons do not modify the leading result 
obtained in the context of the cold plasma theory. The derivation of the dispersion relations 
in the case of a warm plasma can be performed, for instance, within a kinetic approach where 
it can be shown, following the same calculation discussed in the flat space case [12], that the 
first correction to the leading cold plasma calculation can be recast in an effective redefinition 
of the plasma frequency for the electrons, namely u pc — ► 4 ™°^ where 7 = (1 — (f 2 )) -1 / 2 and 
(v 2 ) is the thermal average of the electron velocity. 
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